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The various nesting and pairing instabilities of the gen- 
eralized Van Hove scenario can be classified via an SO (8) 
spectrum-generating algebra. An SO(6) subgroup is an ap- 
proximate symmetry group of the model, having two 6- 
dimensional representations ('superspins'). This group con- 
tains as subgroups both the SO(5) and SO(4) groups found 
by Zhang, while one superspin is a combination of Zhang's 5- 
component superspin with a flux phase instability; the other 
includes a charge density wave instability plus s-wave super- 
conductivity. This is the smallest group which can describe 
both striped phases and superconductivity. 



Two groups play important roles in understanding a 
Hamiltonian: the symmetry group allows a classification 
of its degenerate eigenstates, while the Lie group of the 
spectrum- generating algebroB (SGA) can be used to an- 
alyze the complete spectrum. SGA's have proven to be 
useful in the study of collective modes in nuclear and 
high-energy physics, while in condensed matter physics 
they have been used to study phase transitions in liq- 
uid He and in one-dimensional (ID) metalso. In the ID 
metals, the SGA is SU(8), with 63 elements and 56 pos- 
sible order parameters including superconductivity and 
charge or spin density (CDW/SDW) waves. This alge- 
bra has also been applied to the two-dimensional (2D) 
Hubbard modeB. However, we show that even for the 
generalized Hubbard model appropriate to the general- 
ized Van Hove scenario, the appropriate SGA is SO (8), a 
considerably smaller algebra. This algebra has a natural 
subalgebra, SO (6), which acts as an approximate sym- 
metry group generalizing Zhangta, and including both his 
SO (5) and SO (4) as subgroups. We identify SO (6) as 
the smallest group which is capable of describing striped 
phases as well as superconductivity. n 

For a one-dimensional (ID) metala, nesting involves 
the two points of the Fermi surface, at ifc^, with Uf the 
Fermi momentum. Since this breaks momentum conser- 
vation (the states +kp and —hp become inequivalent), 
the full group SU(8) must be taken as the SGAcl. On 
the other hand, in two dimensions, the dominant nesting 
arises at Q = connecting the two Van Hove sin- 

gularities (VHS's)Q at (it, 0) and (0, n). Since the points 
±(7r, 0) are equivalent points of the reciprocal space lat- 
tice, nesting singularities involve only order parameters 
even in k. Hence the SGA is a proper subgroup of SU(8) 

the SO (8) algebra of Table I. (Note that there is some 
ambiguity in defining a SGA: here we define it as the 
algebra which contains the mean-field Hamiltonian.) 
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FIG. 1. Matrix Representation of SO(8), using the short- 
hand o± = o±oK 

There is a combinatoric interpretation of this SO (8) 
which is independent of any particular hamiltonian. Con- 
sider an electronic system with a two-fold orbital degen- 
eracy (labelled 1, 2) in addition to the spin degeneracy. 
The four creation operators can be written as C\p CL, 
Cj\, and Cvji. Including both particle-hole (C^C) and 

particle-particle (C^C^ or CC) operators, there are 28 
pair operators, whose components define the Lie algebra 
of SO (8) [recall that for SO(N), the Lie algebra contains 
N(N — l)/2 elements]. Particular linear combinations of 
these elements are listed in Table I. Figure [l rewrites 
these elements as an explicit representation of the Lie al- 
gebra of 50(8). The 28 generators are the antisymmet- 



ric matrices L lJ , with matrix elements = 8\8l 



Figure [j] illustrates the equivalences as the lower half of 
an antisymmetric L-matrix. The operators satisfy the 
Lie algebra, with standard 5*0(8) commutation rules 



[L ij ,L k 



8jmL l 



5imL jk - 8jkL lm ). 



(1) 



50(8 — M) subalgebras can be formed by eliminating 
M rows of the L matrices, along with their correspond- 
ing columns. These will be designated as ...,1m}, 
where I\, ...,Im are the rows (and columns) which have 
been eliminated. For .instance, {234} is the £0(5) alge- 
bra studied by Zhanga _ 

In a generalized Hubbard modeD, the creation opera- 
tors can be expanded in terms of operators localized near 
the corresponding VHS's: 



I((-l)-^L(r) + (-ir^L(r)), 



(2) 



where ^>\ a and ip\ a are slowly varying functions of po- 
sition r = a(xi,yi). A more precise definition is given 
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Table I: Generators of SO (8) Lie Algebra 
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Table II: Interaction Terms 

I (G2^G 3 )SrA^(r)A d (r) I 

{2G l +G 3 -G i )^[Ocpw{r)Y 
(G 3 + G4-2G 1 )£HOjc(r)]" 

(G i ^G 3 )^Oj S (r f ) -Ojsj^ 

\-(G 3 + G 4 )S f Og gtv (r) ■ dspw(r)\ 

in Refs.i'i. The Lie algebra of Table I corresponds to 
O -> SfO(r), with C\ a -> Vt(^), etc. With this defini- 
tion, the operators become equivalent to those introduced 
by SchulzQ and ZhangQ. The SGA G is defined in Fourier 
space as G = ® gg^, where gj is the algebra of a particu- 
lar fc-component of the Fourier transformed operators of 
Table I. 

The interaction terms in the generalized Hubbard 
hamiltonianQ can be written in terms of pairs of these 
operators, Table II. Here the Gi's are coupling constants, 
which can be related to the Hubbard U and to various 
near-neighbor interaction termsQ. For the pure Hubbard 
model, G\ = G2 = G 3 = G4 = U /Airt and t is the nearest 
neighbor hopping parameter. The form of the interaction 
term is not unique, since a number of alternative terms 
can arise by anticommuting the operators. 
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FIG. 2. 50(4) weight diagram of V, where S z (J z ) is the 
z component of spin (pseudospin). 

While SO (8) is the SGA of the generalized Hubbard 
model, Zhang's SO (5) is an (approximate) symmetry al- 
gebra of the same model - in particular, the collective 
modes (SDW and d-wave superconductivity for SO (5)) 
are elements of the SGA, but are components of a su- 
perspin, which transforms under the symmetry group. 
For the generalized Hubbard model, the natural (ap- 
proximate) symmetry group is SO (6), defined as follows. 
50(8) contains an 50(3) algebra, which we call isospin, 
generated by T — t,T± — Ocdw±Ojc- This algebra is 
the algebra of the VHS's: the z-component of the isospin, 
To, measures the excess population of the 1-VHS over the 
2-VHS. Those operators which do not commute with To 
can lead to a nesting or pairing instability. Hence, the im- 
portant transformation group of the VHS's is the 50(6) 
subgroup {23} which commutes with To, leading to the 
decomposition scheme: 

SO(8)^SO(6) {23} eV + ©V_©r, (3) 

under which 28 -> (15, 0) + (6, l) + (6, -1) + (1,0), where 
(m, n) denotes representation m of 50(6) and eigenvalue 
n of T . The 6- vectors can be denoted V± = V± V with 

V = {L21, L42, Lq2, L72, -^82} 

V = {L 3 i, L43, L 53 , L 6 3, £737 ^83}, (4) 

shown boxed in Fig. |l|. The group structure of V is 
shown in Fig. where J z is the z-component of the 
pseudospin operator introduced by Yang and ZhangtHI 
and S dw± = T(Osdw^ ± iOsDWy)/V2- An analo- 
gous diagram can be drawn for V. The group 50(6){23} 
transforms the components of each of these 6-vectors 
among themselves, without mixing the two vectors, while 
r transforms the vectors into each other. 
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A number of points should be noted. (1) The SO(6) 
group {23} contains Zhang's SO (5) group as a subgroup, 
as well-as the the 50(4) group introduced by Yang and 
Zhango. Moreover, V, Fig. |[ combines Zhang's SO (5) 
superspin with_0 jc, which is essentially equivalent to 
the flux phaseEi 

(2) The twelve components of superspin are precisely 
the collective modes identified earlier by SchulzOI, and 
most of them have been found to play an important role 
in the cuprates: s-wave superconductivity in electron- 
doped and (possibly) overdoped cuprates, CDW's-Jiear 
optimal dopingiI30, the flux phase near half fillingli-S't^l. 

(3) For a bare band dispersion (neglecting interactions) 
of the form 

- 4t cos k x a cos k y a, (5) 



— 2i(cos k x a + cos k y a) - 



two parameters control the symmetry of the quadratic 
part of the generalized Hubbard hamiltonian, t' and 
jl = Ep — Ey, the shift of the Fermi level Ep from the 
VHS Ev ■ When both parameters are zero (half filling 
with square Fermi suriace) the hamiltonian has an extra 
pseudospin symmetrylMx In this case, the nature of the 
ground state instability is controlled solely by the inter- 
action terms (the G's). A pure Hubbard interaction (U) 
breaks the SO(6) symmetry (Table II): 



SO(6) -> SO(3)®SO(3) 



(6) 



with one -SjO(3) ordinary spin, and the other the 
pseudospinEI Both 6-vectors are broken down to pairs 
of 3- vectors 



V^{0 JSi }®{At,A s ,0 CDW } 
V -> {OsDWi} © {A d , A d , Ojc}; 



(7) 



however, there remains an accidental degeneracy of one 
vector (Ojs) with the opposite pseudovector. At half 
filling the lowest energy state is Osdw- As discussed 
below, this weak coupling result must be corrected for 
strong correlation effects. 

When the Fermi surface is distorted from square, either 
by doping away from half filling (/}) or by introducing 
second-neighbor hopping terms t', the pseudospin degen- 
eracy is broken, in such a way as to favor pairing onec 
nesting instabilities. ThiSfCanj.be seen by Hartree-Fockt3 
or renormalization groupuOE3|-apalyses or by a linear 
response analysis (following RefJla). 

If the superspin is written as O (a twelve component 
vector incorporating both representations), then in lin- 
ear response theory it is assumed that there is an ap- 
plied field ho (also a 12- vector) which couples to O. The 
hamiltonian in the presence of ho is 



H = ^ioH a \ a Ha + ho-0 



(8) 



(here the terms in G have been neglected) , with resulting 
free energy 



F (O) = n ( t i,O,T) + l iN 
fi (/x, O, T) = -2k B TE^ ln(l + e^*-^'^) 

+ho ■ O, (9) 

with Er the quasiparticle energy found by applying a 
Bogoliubov.J/alentin transformation to Eq. || or via SGA 
techniquesErD. The expectation value of each superspin 
component Oi is found from 



dh~ 



Oi 



0. 



and the corresponding susceptibility is 
Xch = Jim (t-*-). 

h o ^0 n Oi 



(10) 



(11) 



Including the interaction terms, the Hartree-Fock free 
energy becomes 



O 



(12) 



leading to an instability of the ith mode when 

1 + 2xoiGi = 0. (13) 

If the quadratic hamiltonian is symmetric under SO(6), 
then the component with the most negative Gi is the 
first to diverge. For finite jl or i',the hamiltonian still 
preserves particle number, so there are only two inde- 
pendent susceptibilities, the particlc-hole susceptibility 
Xoo and the pair susceptibility X02, with 



Xoo 



-2Se 



e ko e k+Q,a 



X02 



ka 



(14) 



Note that in nearest-neighbor hopping models [t 1 = 0) 
-eg , and the two expressions become equiva- 



C k+Q,<r 



lent when ep = - i.e., at half filling. 

Figure |H illustrates the doping dependence of these sus- 
ceptibilities for a Hubbard band with nearest-neighbor 
hopping only (t' = 0) . The point of maximum instability 
(largest x) coincides with the point at which the VHS 




FIG. 3. Susceptibilities xoo (dotdashed line) and X02 (solid 
line) vs. band filling 1 + x for Eq. H with t' = 0. 
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crosses the Fermi level - half filling when t' = 0. When 
x = t' = 0, the susceptibilities are degenerate, xoo = X02, 
as expected from the pseudospin syrrrrnetryll9. However, 
as soon as the system is doped away from half filling 
(x 7^ 0) the electron-hole susceptibility drops precipi- 
tously, whereas the pair susceptibility falls off much more 
gradually. A similar effect arises if the system is main- 
tained at optimal doping (the VHS), but the parameter 
t' is varied - indeed X02 actually increases with increasing 
t' , Fig ||. This striking difference is readily understood: 
the electron-hole susceptibility involves inter- VHS nest- 
ing, which gets progressively worse as the Fermi surface 
gets more curved, whereas the electron-electron suscepti- 
bility involves intra- VHS scattering, and increases with t' 
as the Fermi surfaces become nearly ID near the VHS's. 
(In Figs. [||4], the logarithmic divergence at the VHS was 
cut off by adding a small imaginary term to the denomi- 
nator of x-) 



are associated with a CDW: the low-temperature tetrago- 
nal phase is nearly coterminus with the long-range SDW- 
ordered phase, and the fact that the charged stripes are 
best seen by neutron diffraction suggests a strong associ- 
ated lattice distortion. There is considerable additional 
evidence that phonons and structural instabilities play an 
important role in the doped materialta Hence, for a de- 
tailed description of the doping dependence of the pseu- 
dogap, striped phases, and extended VHS's, it may be 
necessary to recognize that strong electron-phonon cou- 
pling can lead rto a crossover to a groundstate involving 
the V 6-vectonl3. 

MTV's research is supported by the Dept. of Energy 
under Grant # DE-FG02-85ER40233. Publication 722 
of the Barnctt Institute. 




FIG. 4. As in Fig. |, but with 4t' = E F . 

Figure His consistent with the RG results of SchulzB. 
From Eq.TL3L when all x' s are equal the order parameter 
associated with the most negative G is the first to go sin- 
gular. For the pure Hubbard model, this means that the 
leading instability at half filling is the SDW. When the 
material is doped, the SDW susceptibility plummets, and 
at some point d-wave superconductivity becomes favor- 
able. In agreement with Zhangp, the shift of the Fermi 
energy from the VHS is a relevant parameter in driving 
this SDW^ d-wave superconducting transition. 

Despite the simplicity of this picture, a purely SO (5) 
model cannot explain the full physics of the cuprates. 
First, the above analysis is in the weak coupling limit, 
and a strong coupling reanalysis of Table II (J-term 
dominant) shows that the flux phase - not included in 
SO (5) - is the lowest energy state. A second problem is 
in La2- a; Ba a ;Cu04 and L^-z-yNdySr^CuO^ near x = 
1/8, where the striped phase is commensurately pinned, 
leading to long-ranged magnetic and charge orderEj. At 
the same time, superconductivity is strongly suppressed, 
demonstrating that whatever the driving force for charge 
order may be, it is not superconductivity, but is in com- 
petition with superconductivity. Since SO (5) only allows 
for antiferromagnetism and superconductivity it docs 
not have sufficient flexibility to properly describe this sit- 
uation. There are strong hints that the charged stripes 



1 "Dynamical Groups and Spectrum Generating Algebras", 
ed. by A. Bohm, Y. Ne'eman, and A.O. Barut (Singapore, 
World Scientific, 1988). 

2 A.I. Solomon and J.L. Birman, J. Math. Phys. 28, 1526 
(1987). 

3 A.B. Eriksson, T. Einarsson, and S. Ostlund, Phys. Rev. 
B52, 3662 (1995). 

4 S.-C. Zhang, Science 275, 1089 (1997). 

5 J. Solyom, Adv. Phys. 28, 201 (1979). 

6 T.M. Rice and G.K. Scott, Phys. Rev. Lett. 35, 120 (1975). 

7 H.J. Schulz, Phys. Rev. B39, 2940 (1989). 
8 R.S. Markiewicz, Physica C210, 235 (1993). 

9 R.S. Markiewicz and M.T. Vaughn, to be published, J. 
Phys. Chem. Sol. (cond-mat/9709137), and unpublished. 

10 C.N. Yang and S.-C. Zhang, Mod. Phys. Lett. B4, 759 
(1990); S.-C. Zhang, Phys. Rev. Lett. 65, 120 (1990). 

11 I. Affleck and J.B.Marston, Phys. Rev. B37, 3774 (1988). 

12 R.S. Markiewicz, J. Phys. Chem. Sol. 58, 1179 (1997). 

13 R.S. Markiewicz, Phys. Rev. B56, 9091 (1997). 

14 T.C. Hsu, Phys. Rev. B41, 11379 (1990); R.B. Laughlin, 
J. Phys. Chem. Sol. 56, 1627 (1995); X.-G. Wen and P. A. 
Lee, Phys. Rev. Lett. 76, 503 (1996). 

15 C.A. Balseiro and L.M. Falicov, Phys. Rev. B20, 4457 
(1979). 

16 I. E. Dzyaloshinskii, Zh. Eksp. Teor. Fiz. 93, 1487 (1987) 
[Sov. Phys.-JETP 66, 848 (1987)]. 

17 H. Schulz, Europhys. Lett. 4, 609 (1987). 

18 T. Moriya, "Spin Fluctuations in Itinerant Electron Mag- 
netism" (Springer, Berlin, 1985). 

19 J.M. Tranquada, B.J. Sternlieb, J.D. Axe, Y. Nakamura, 
and S. Uchida, Nature 375, 561 (1995); J.M. Tranquada, 
J.D. Axe, N. Ichikawa, A.R. Moodenbaugh, Y. Nakamura, 
and S. Uchida, Phys. Rev. Lett 78, 338 (1997). 



4 



